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Abstract. We study power density spectra (PDS) of X-ray flux variability in binary systems where the accretion 
flow is truncated by the magnetosphere. PDS of accreting X-ray pulsars where the neutron star is close to the 
corotation with the accretion disk at the magnetospheric boundary, have a distinct break/cutoff at the neutron star 
spin frequency. This break can naturally be explained in the "perturbation propagation" model, which assumes 
that at any given radius in the accretion disk stochastic perturbations are introduced to the flow with frequencies 
characteristic for this radius. These perturbations are then advected to the region of main energy release leading to 
a self-similar variability of X-ray flux P oc / — 1 - - -1 - 5 . The break in the PDS is then a natural manifestation of the 
transition from the disk to magnetospheric flow at the frequency characteristic for the accretion disk truncation 
radius (magnetospheric radius). The proximity of the PDS break frequency to the spin frequency in corotating 
pulsars strongly suggests that the typical variability time scale in accretion disks is close to the Keplerian one. 
In transient accreting X-ray pulsars characterized by large variations of the mass accretion rate during outbursts, 
the PDS break frequency follows the variations of the X-ray flux, reflecting the change of the magnetosphere size 
with the accretion rate. Above the break frequency the PDS steepens to ~ f~ 2 law which holds over a broad 
frequency range. These results suggest that strong / aperiodic variability which is ubiquitous in accretion 

disks is not characteristic for magnetospheric flows. 
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1. Introduction 



It has been recognized since the beginning of X- 
ray astronomy that the flux of accreting X-ray bina- 
ries demonstrates strong aper i odic variability ( s ee e.g . 
Rappaport. Doxsev. fe Zaumenl ll9Tlh lOda et aT 1974 ). 
Almost immediately after the discovery, the noise in the 
X-ray light curves of accreting binaries (like e.g. Cyg X- 
1) was explained as a superposition of randomly occur- 
ring X-ray emissio n flashes (sho ts) of similar duration (the 
shot noise model, This provided an expla- 

nation to the shape of the power density spectra (PDS, 
the Fourier transform of the autocorrelation function of 
the lightcurve of a source) of different X-ray sources. 

However, the accumulation of more data posed se- 
rious questions to this paradigm. In particular, it was 
very hard to explain a huge range of the X-ray vari- 
ability time scales observed in some s ources (see e.g. 
Churazov. Gilfanov. fc Revnivtsevl 120011 ) and the linear 
correlation between th e variability amplitude an d the av- 
erage flux of sources (jUttlev fc McHardvl 120011 ). Indeed, 
to explain observed large variability amplitude, the in- 



dividual flashes/shots in the shot noise model should be 
very powerful. Therefore these flashes must come from 
the innermost region of the accretion flow, where most 
of the energy is released. The characteristic time scales 
in this region are very short - milliseconds or tens of 
milliseconds for stellar-mass compact objects. However, 
very often, e.g. in the soft/high state of accreting bina- 
ries, the observed power spectra have a power law shape 
extending down to frequencies as low as 10 -5 — 10~ 6 Hz, 
i.e. 5-7 orders of magnitude longer time scales than all 
time scales char acteristic for the region of the main en 
ergy release (seelChurazov. Gilfanov. fc Revnivtsevll2001 



Gilfanov fc Arefievl 120051 ). 

A very promising 
X-ray variability of 
"pert urbat ion propagation" model (iLvub arskii 



model for the 
accreting sources 



aperiodic 
is the 



1997 



Churazov, Gilfanov, & Revnivtsev 



2001 



Kotov. Churazov. fc Gilfanovl 2001 
20061 ). In this model, the X-ray flux variability is caused 



Arevalo fc Uttlev 



by the variations of the instantaneous value of the mass 
accretion rate in the inner accretion flow. In turn, the 
variations of the mass accretion rate are due to the 
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perturbations introduced to the accretion flow by the 
stochastic variations of the disk viscous stresses. In 
this model the observed variability is a multiplicative 
superposition of perturbations introduced at different 
radii. Assuming that the fractional amplitudes of the 
mass accretion rate perturbations are the same at all 
radii, the PDS of the emerging lightcurve will nat- 
urally appear as a self-similar power-law with slope 
— 1... — 1.5 up to the maximal frequen cies that can be 
generated in the disk (|Lvubarskiilll997h . Direct magneto 
hydrodynamic sim u lation s of accretion flows ( se e.; 



Brandenburg et al. 19951: Balbus fc Papaloizou 199 



Hirose. Krolik. fc Stonei 2006 ) provide further support to 
this semi-phenomenological model. In particular, these 
simulations show that perturbations in the instantaneous 
mass accretion rate generated at any given radius of the 
disk have characteristic time scale proportional to the 
local dynamical time. 

This model of the aperiodic X-ray flux variability im- 
plies that the presence of the accretion disk edges, both 
outer and inner, should be reflected in the noise proper- 
ties of the X-ray light curve. S ignatures of outer edges o f 
accretion disks were found bv iGilfanov fc Arefiev ( 20051 ) 
in the low frequency parts of the noise power spectra of 
low mass X-ray binaries. 

Accretion disks around compact stars in X-ray bi- 
naries should also have inner edges, which should 
manifest itself in the power spectrum of their X-ray 
light curves. Specifically, a definite break is expected 
to be present in the PDS at the characteristic fre- 
quency of variability generated at the inner edge of 
the disk. At frequencies below this break the power 
spectrum is expected to be produced in the accre- 
tion disk and have a self-similar slope about -1.0...- 
1.5 (Lvubarskiil Il997t Churazov, Gilfanov. fc Revnivtsev 
200lt IGilfanov fc Arefievl 120051 : lllevnivtsev et all 120061 ) 
while at higher frequencies the character of the flow 
changes and the PDS slope may be different. 

In accreting X-ray pulsars and intermediate polars the 
central compact object (a neutron star or a white dwarf) 
has a strong magnetic field which can disrupt the disk- like 
accretion flow at the magnetospheric boundary, or even 
prevent the formation of accretion disk at all (like in po- 
lars), dividing the flow into two distinct parts - the accre- 
tion disk and the magnetospheric flow. The noise proper- 
ties of these flows may be very different. 

In the present paper we compare PDS of different types 
of accreting X-ray binaries and discuss the observational 
support for the qualitative picture of aperiodic variability 
outlined above. 



— X-ray binaries with compact objects, which have mag- 
netospheres powerful enough to disrupt the accretion 
disk at large distances (accreting X-ray pulsars, inter- 
mediate polars) 

— X-ray binaries with large magnetospheres in which the 
accretion disk does not form and the accretion pro- 
ceeds along the magnetosphere from the very begin- 
ning (polars) 

In a fair fraction of persistent accreting X-ray pulsars 
the spin period of neutron star is observed to be close 
to synchronization (corotation) with the Keplerian rota- 
tion of the accretion disk at the magnetospheric bound- 
ary, which is explained by the standard description of the 
interaction of accretion disk with a magnetized neutron 



star (see e.g. iDavidson & Ostrikerl 


19731: IShakura 1975: 


Litmnov & Shakura! 


197a iGhosh & Lambl Il979t Corbetl 


1984: Ziolkowski 


1985). 



On the contrary, due to much larger moments of iner- 
tia of white dwarfs, it takes a much longer time to bring 
accreting magnetic white dwarfs in cataclysmic variables 
into corotation with surrounding accretion disks, and their 
magnetospheres rotate with periods much larger than the 
Keplerian ones at the inner disk edge. 

The perturbation propagation model makes distinct 
predictions for the noise properties of accretion flows in 
these classes of objects. Namely: 

1. The noise power spectra of accreting sources with large 
magnetospheres and without accretion disk (polars, 
e.g. AM Her), in which matter is transferred directly 
from the companion star to the compact star via mag- 
netospheric accretion, should be different from those 
of sources with accretion disks and without magneto- 
spheres (e.g. accreting black holes like Cyg X-l in their 
soft state). 

2. The noise power spectra of sources with accretion 
disk truncated at the magnetospheric boundary should 
have a break corresponding to the characteristic fre- 
quency in the disk near its inner boundary. 

3. If the size of the magnetosphere changes (e.g. as a re- 
sponse to the change in mass accretion rate like in 
transient pulsars), the break frequency must change 
correspondingly. 

4. For persistent X-ray pulsars which are close to coro- 
tation, the comparison of the PDS break frequency 
with the compact object spin frequency can be used 
to determine the ratio of the characteristic frequency 
of perturbations generated in the disk to the Keplerian 
frequency. 



2. Truncated accretion disks in different classes of 
sources 

We shall consider several classes of accretion X-ray bina- 
ries: 



2.1. Breaks in the power spectra spectra of 
magnetized accretors 



It has been 



Hoshino fc Takeshima 



noticed in earlier studies (e.g. 

1993) that noise power spec- 



tra of accreting X-ray pulsars typically have breaks 
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Fig. 1. Power density spectra of persistent accretion pow- 
ered X-ray pulsars. The frequency shown along the X- 
axis is expressed in units of the compact object spin fre- 
quency. The power density spectra are multiplied by fre- 
quency to show the square of the fractional RMS per 
decade of frequency. For comparison with the observed 
profiles, the thick dashed curve shows an analytical model 

p oc r 1 (i+ (///o) 2 )- - 5 . 



around the pulse frequency. This is shown more clearly 
in Fig[TJ where we plot the noise power spectra of several 
persistent X-ray pulsars, including Cen X-3, 4U 1626-67, 
Vela X-1, 4U 1538-52, GX 301-2, X Persei, GX 1+4, 
and of the magnetic white dwarf (polar) AM Her. The 
frequency scale of their power spectra is normalized 
through multiplication by the pulse period of the sources. 
The variability, associated with regular pulsation was re- 
moved from the original X-ray light curves by subtracting 
folded segments of light curves with a duration of 10-20 
spin periods . Strictly speaking this procedure does not 
remove the contribution of regular pulsations completely 
since often the periodic and aperiodic variabilities show 
signs of nonlinear inte ractions (see e.g. the discussion in 
Tsygankov et al.ll2007l ): howe ver, it is good enough as the 
first approximation (see e.g. iFinger. Wilson, fc Harmonl 
19961) . 

It is seen that all power spectra of accreting X-ray 
pulsars show clear breaks approximately at their spin fre- 
quency. PDSs of all sources have a similar power-law slope 
above the break frequency irrespective of the PDS form 
below it. 

In the framework of the pertu r bation 
propagation model 



truncation of the accretion disk flow and its conversion 
into a magnetospheric flow. 

The proximity of the PDS break frequency to the spin 
frequency (which in the case of corotating systems is close 
to the frequency of the Keplerian rotation around the com- 
pact object at the inner edge of the accretion disk) allows 
us to conclude that the characteristic time scale of vari- 
ability produced at some distance from the central com- 
pact objects are close to the local Keplerian time scale. 

Fig.[T]shows that above the break frequency the power 
spectrum of the flux variability typically is a power law 
with the slope close to ~ —2. We can not exclude that 
such slope of the flux variability is a property of the mag- 
netospheric accretion, however, we note that similar slopes 
of power spectra sometimes observed in accreting systems 
without magnetospheric accretion, for example unmagne- 
tized cataclysmic variables (see e .g. Kraicheva et al. 19991 
Pandel. Cordova, fc Howelill2003l ) . 



sec 



c. 



Churazov. Gilfanov. fc Revnivtsevl I2001I ). we can in 



Lvubarskiil 11997 



terpret this observational fact as a signature of the 



2.2. Change of the magnetospheric size with mass 
accretion rate 

If the break frequency in the noise power spectra of ac- 
creting X-ray pulsars indeed reflects the time scale of the 
noise generation at the inner boundary of the accretion 
disk/flow, its value should depend on the mass accretion 
rate in the binary system. Increase in the mass accretion 
rate decreases the size of the magnetosphere (and hence 
the inner radius of the disk) and brings the system off 
the corotation, so the characteristic frequency at the inner 
edge of the disk/flow increase. A similar situation should 
be observed in the case of luminous intermediate polars 
(accreting magnetized white dwarfs with moderate mag- 
netospheres) which typically are always out of the corota- 
tion. 

We can verify this hypothesis by examining power 
spectra of bright transient X-ray pulsars (like A0535+26, 
41X0115+63, V0332+53, KS 1947+300), which demon- 
strate wide range of X-ray luminosities during the out- 
bursts, and power spectra of luminous intermediate polars 
(e.g. V1223 Sgr). 

In Fig|5] we show examples of power spectra of the 
pulsar A0535+26 during its bright outburst in 2005 and 
the power spectrum of V1223 Sgr. The power spectra of 
A0535+26 were averaged over two time intervals: a) with 
small X-ray luminosity (L x ~ 10 36 erg/s) (low mass ac- 
cretion rate) and b) with higher luminosity (L x ~ 10 37 
erg/s) (higher mass accretion rate). The power spectrum 
of V1223 Sgr is averaged over all publicly available RXTE 
observations. 

From Fig [2] it is clear that the power spectra of 
A0535+26 in different luminosity states differ significantly 
at the frequencies higher than the pulse frequency, while 
at lower frequencies the PDSs are almost identical. The 
difference at high frequencies is essentially an addition of 
an extra noise component, which is presumably generated 
in the ring of accretion disk between the radii correspond- 
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ing to the magnetosphere size in the high accretion rate 
(small radius) and low accretion rate (large radius). This 
ring and associated variability were absent in the state 
with low accretion rate. 

The noise power spectrum of V1223 Sgr closely resem- 
bles that of A0535+26 in the bright state. This similar- 
ity reflects the fact that both systems are out of corota- 
tion with their accretion disks - the magnetospheres are 
squeezed by the increased accretion flow, and the inner 
parts of the accretion disks rotate much faster than the 
central object. 



0.01 r 
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Fig. 2. Power density spectra of accreting X-ray pulsar 
A0535+26 at low accretion rate (labeled as "faint") and 
in the strong spin-up regime during the outburst (labeled 
as "bright"). The upper curves shows the power spectrum 
of luminous intermediate polar V1223 Sgr 
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Fig. 3. Dependence of the break frequency in the noise 
power spectrum of A0535+26 on the 3-20 keV X-ray 
flux (filled circles). The dependence of the QPO fre - 
quency observed by Finger. Wilson, fc Harmonl ( 19961 ) 
during the source outburst in 1994 is shown by open 
dotted circles. The solid circles show the results of 

when the QPO fre- 



Finger. Wilson, fc Harmonl (|199(: 
quency is recalculated to the break frequency as fbreak = 
2.5 x Jqpo (see Fig. [U the upper panel). The dashed line 
shows the prediction of the simplest "magnetospheric" 
model of the break frequency /b oc £ x 3 ^ 7 described in 
the text. The dotted line shows the neutron star spin fre- 
quency. 



If the characteristic frequency /o of the noise at the 
magnetospheric boundary i? m is proportional to the fre- 
quency of the Keplerian rotation j/k of matter at the in- 
ner edge of the accretion disk i?; n w i? m , we can relate 
the observed break frequency to the instantaneous value 



of the mass accretion rate M (see e.g.lPrin glcfc Rec: 

~3 



Lamb. Pethick. fc Pinesl 1973t Davidson fc Ostrike 



1972; 



1973; 



Bildsten et al.lll997l ): 

27WK 



{GMY' 2 R^' 2 



where 



i?m » M 4/7 (GM)-V7M-2/7 



is the standard expression for the magnetospheric radius, 
fj, is the dipole magnetic moment of the neutron star. 
Therefore we can anticipate that the break frequency will 
follow the dependence: 



(1) 



It is exactly what we see during the evolution of the 
outburst of A0535+262 observed by RXTE in 2005 (we 



have substituted the X-ray luminosity L x 
stead of the mass accretion rate M into Eq. ([1])). The 
dependence of the break frequency (ft, in Eq. (p}) on the 
X-ray flux of A0535+262 is shown in Fig. [3] (filled circles). 
For this this plot we have used all observations of RXTE of 
this outburst after ~MJD 53613, when the stable regime 



of accretion was es tablished (see ICaballero et al.l 12008 



Postnov et all 2008) 



3/7 

It is interesting to note that the dependence / oc L x 
was previously established for the centroid frequency of 
quasi-periodic X-ray oscillations in the power spectra of 
A0535+26 detected during its g iant outburst in 1994 
( Finger. Wilson, fc Harmonl Il996j ). The authors argued 



that the oscillations originate at the inner boundary of 
the accretion disk and are related either to the Keplerian 
rotation at the inner edge of the disk or to the beat fre- 
quency between the Keplerian r otation and that of th e 
neutron star magnetosphere (e.g. Alpar fc Shahaml 1985 ). 

Actually we can try to combine our present 
RXTE measurements with QPO studies by 
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Fig. 4. Upper panel: power spectra of Cen X-3 and 
4U1626-67 with low-frequency QPOs (the variability as- 
sociated with regular pulsations has been subtracted from 
the original X-ray light curves). The frequency axis has 
been scaled by the break frequency in the corresponding 
power spectrum. Lower panel: the broadband X-ray en- 
ergy spectrum of A0535+26 observed by RXTE during 
its outburst in 2005. For this spectrum the ratio of en- 
ergy fluxes in 20-100 keV and 3-20 keV energy band is 
approximately ~0.55. 



Finger, Wilson, fe Harmon! (|1996 
facts: 



making use of two 



a) when low-frequency QPOs are observed in power spec- 
tra of accreting pulsars, their centroid frequency is re- 
lated to the br eak frequency as ~1:2.5 (F i gHl u pper 
panel, see also lAngelini. Stella, fc Parmail p989) for 
QPO in EXO 2030+375). The correlation of the QPO 
and the break frequencies is also observed in LMXBs, 
although in these objects the centroid of the QPO 
feature is usually above the br eak frequency (see e.g. 
Wiinands fc van der Kliill999h . 



b) during RXTE observations of the source outburst in 
2005, the X-ray flux of the A0535+26 in the energy 
band 20-100 keV (the range wh ere measurements by 
Finger. Wilson, fc Harmo nl ll996l were done) is approx- 



imately factor of 0.55 of the X-ray flux in the energy 
band 3-20 keV (see the RXTE spectrum shown in the 
lower panel of FigU]) . 

The dependence of the QPO frequency on the X-ray flux 
during its 1994 outburst, renormalized using the above 
factors, is shown in Fig. [3] (solid open circles). The renor- 
malized dependence (solid open circles) perfectly contin- 
ues the observed break frequency-flux dependence (filled 
circles) and corresponds to Eq. ([TJ) (the dashed line). 

Prominent QPO features simila r to those detected 
by iFinger. Wilson, fc Harmonl (|l996l ) are not always ob- 
served in the power spectra of accreting X-ray pulsars. 
On the other hand a break in the PDS is more ubiquitous 
and therefore the diagnostics of the accretion flow based 
on the break frequency can be applied to larger datasets. 
For example, the break frequency in the noise power spec- 
trum can be used as an estimate of the dipole magnetic 
moment of compact objects using Eq. (fTj). 

3. Conclusions 

We studied aperiodic variability of the X-ray flux from 
accreting binaries, in which the truncation of the disk-like 
accretion flow by the magnetosphere of the compact object 
is important. The results can be summarized as follows: 

— There is a distinct break in the Power Density Spectra 
of accreting magnetized neutron stars and white 
dwarfs, apparently associated with the change of the 
disk-like accretion flow to the magnetospheric flow 
near the Alfvenic surface. 

— In transient systems with variable X-ray luminosity the 
PDS break frequency fb changes with the X-ray lumi- 
nosity (mass accretion rate) as fi, rx L^ 7 , in agreement 
with the standard theory of accretion onto magnetized 
compact stars. 

— This break can naturally be explained in the "per- 
turbation propagation" model, which assumes that at 
any given radius in the accretion disk stochastic per- 
turbations at frequencies characteristic for this radius 
are introduced to the flow. These perturbations are 
then advected by the flow to the region of main en- 
ergy release leading to a self-similar form of the PDS 
P rx The break in the PDS corresponds to 
a frequency characteristic for the accretion disk trun- 
cation radius (the magnetospheric radius). 

— We suggest that the PDS break frequency is directly 
related to the magnetospheric radius for a given value 
of the mass accretion rate and can be used to estimate 
the magnetic moment of accreting compact stars. 

— For systems which are close to corotation (accreting 
X-ray pulsars) the PDS break frequency is close to the 
spin frequency of the neutron star. This strongly sug- 
gests that the characteristic frequency of perturbations 
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introduced to the accretion flow in the disk is of order 
of the local Keplerian frequency. 
— In all studied objects the PDS above the break fre- 
quency follows the P ~ f~ 2 law over a broad range of 
frequencies, suggesting that strong f- 1 —- 1 - 5 aperiodic 
variability which is ubiquitous in accretion disks is not 
characteristic for magnetospheric flows. 
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